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ABSTRACT 

We have determined the ratios of total to selective extinction directly from observations in 
the optical V band and near-infrared J band toward the Galactic center. The OGLE (Optical 
Gravitational Lensing Experiment) Galactic bulge fields have been observed with the SIRIUS 
camera on the IRSF telescope, and we obtain Ay/Ey^j = 1.251 ± 0.014 and Aj/Ey-j — 
0.225 ±0.007. From these ratios, we have derived Aj/A v = 0.188 ±0.005; if we combine Aj/Ay 
with the near-infrared extinction ratios obtained by Nishiyama et al. for more reddened fields 
near the Galactic center, we get Ay : Aj : Ah ■ Ak s — 1 : 0.188 : 0.108 : 0.062, which implies 
steeply declining extinction toward the longer wavelengths. In particular, it is striking that the 
K$ band extinction is ~ 1/16 of the visual extinction Ay, much smaller than one tenth of Ay 
so far employed. 

Subject headings: dust, extinction — stars: horizontal-branch — Galaxy: center 



1. INTRODUCTION 

The wavelength dependence of interstellar ex- 
tinction provides important diagnostic informa- 
tion about the dust grain properties. Interstel- 
lar extinction law shows a large range of variabil- 
ity from one line of sight to another, especially in 
the ultraviolet and optical wavelengths. In com- 
paring the wavelength dependence among different 
lines of sight, the normalization of the extinction 
curves by the total extinction Ay, instead of the 
usual color excess Eb-v, is vitally important, as 
Cardelli, Clayton, & Mathis (1989, CCM) have 
shown. According to CCM, the variation in Galac- 
tic extinction curves from the ultraviolet to the 
optical is described by a single parameter, which 
itself is the ratio of total to selective extinction 
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Rv = Av/Eb-v- Thus the ratio Ry, or more 
generally R\ — A\/ E\>_\, is very important, al- 
though very difficult to obtain. An usual way to 
determine Ry is to extrapolate the ratio of color 
excesses E\_y / Eb-v to A = oo with reference to 
some model, and can be compromised by emission 
or scattering by dust grains near the stars. 

The past decade has seen a new method to de- 
termine the ratio of total to selective extinction 
R\. The method was first proposed by Wozniak 
& Stanek (1996); in essence, it simply measures 
the regression of the mean color of red clump (RC) 
stars in the Galactic bulge on their mean magni- 
tude (RC method). The increase in the amount 
of dust in a line of sight causes the clump fainter 
and redder. The slope of these changes in a color- 
magnitude diagram (CMD) is equivalent to R\. 
The method has been developed by Stanek (1996), 
Udalski (2003), and Sumi (2004) in the V and / 
bands 1 , and recently applied to the near- infrared 
wavebands by Nishiyama et al. (2006a, hereafter 
paper I). 



They regard their / filter as similar to the Landolt Cousins 
J, whose effective wavelength is about O.SO^m. 
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As a result, a somewhat surprising suggestion 
was first made by Popowski (2000) for the in- 
terstellar extinction toward the Galactic bulge; 
the ratio of total to selective extinction Ryi = 
Av/Ev-i is approximately 2.0, much smaller 
than the ratio 2.5 for the standard (Ry = 3.1) 
CCM extinction curve (Udalski 2003). This was 
confirmed by Sumi (2004) in a larger number 
of fields. Thus the V (0.55^m) to I (0.80yum) 
part of the extinction curve toward the Galac- 
tic bulge seems to be characterized by a smaller 
Rvi- In the near infrared, the extinction law is 
frequently referred to as "universal" , and in fact 
CCM gives an i?y-indcpcndcnt curve. However, 
paper I have shown that the near-infrared extinc- 
tion curve toward the Galactic Center (GC) is 
different from those in the literature, having also 
smaller Ah/Ej_h and Ak s / 'Eh-k s - 

In this paper, we apply the RC method and ex- 
tend the R\ determination to the J band (1.25/um) 
by measuring the J magnitude of the stars whose 
V photometry was obtained in Udalski (2003). 
Furthermore, assuming that the extinction curve 
can be extended to more heavily reddened regions 
observed in paper I, we estimate A\jAy at J, H, 
and Ks wavelengths. 

2. Observation, Data Reduction, and Anal- 
ysis 

2.1. Near-infrared Observations 

All observations in the near-infrared wave- 
lengths were made with the SIRIUS camera (Si- 
multaneous InfraRed Imager for Unbiased Survey; 
Nagashima et al. 1999; Nagayama et al. 2002) at- 
tached to the 1.4 m telescope IRSF (InfraRed Sur- 
vey Facility), on the nights of 2004 May 18 and 
19. The SIRIUS camera provides photometric 
images of a 7'.7 x 7.'7 area of sky in three near- 
infrared wavebands J (1.25/xm), H (1.63^m), and 
Ks (2.14/im) simultaneously. The detectors are 
three 1024 x 1024 HgCdTe arrays, with a scale of 
0'.'45 pixel" 1 . The IRSF/SIRIUS system is similar 
to the MKO system (Tokunaga et al. 2002). 

Over the range of 17 h 52 m < R.A.(J2000.0) < 
17 h 56 m and -30?2 < Dec.(J2000.0) < -29?5, 32 
images in each band were obtained (Fig. 1). These 
fields correspond to the region "A" (( w 0°,i) ~ 
-2°) of Udalski (2003) in the Optical Gravita- 
tional Lensing Experiment (OGLE) II maps of the 
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Fig. 1. — Observed area and fields used for data 
analysis. Thick-line rectangles, labeled as "SC39", 
"SC4" , "SC3" , and "SC37" show the areas ob- 
served by the OGLE project. The fields observed 
by IRSF/SIRIUS are shown by hatched squares. 



Galactic bulge in the V and / bands 2 with an 
absolute photometric accuracy of 0.01-0.02 mag 
(Udalski et al. 2002). The data files containing 
the photometry were downloaded from the OGLE 
website 3 . The Udalski (2003) "A" region consists 
of the five fields SC3, SC4, SC5, SC37, and SC39 
(Fig. 1), but the SC5 data were excluded from 
the following analysis because it suffers from the 
heaviest extinction and its V magnitudes become 
unreliable, as Udalski (2003) also pointed out. 

The observing weather was photometric, with 
seeing of in the J band. Flat fields were 

obtained during each clear evening and morning 
twilight. Dark frames were taken at the end of 
each observing night. A single image comprises 10 



2 In this paper, we do not use the / band data, because the 
OGLE I band filter has a wider long-wavelength wing of 
transmission contrary to the sharper drop of the standard / 
band filter, and the relation between V and / extinction has 
been already studied by Udalski (2003), and Sumi (2004). 

3 http://bulge. princeton.edu/ ~oglc/oglc2/bulgc_maps. html 
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dithered 5 sec exposures. 

The SIRIUS camera provides three (J, H, and 
Ks) images simultaneously, and thus we have the 
data sets of the same region in the three bands. 
However, the uncertainties in the H and Ks pho- 
tometry are larger, and the range of extinction is 
smaller, about a half (H) and a third (Ks) of that 
in the J band 4 . Therefore we have obtained reli- 
able results only in the J band. 
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Fig. 2. — Histograms of the positional differ- 
ence between the OGLE and IRSF coordinates 
for cross-identified stars in R.A. (left) and Dec. 
(right). The values of a obtained by fitting with 
a Gauss function are ft' 088 in R.A. and ft'081 in 
Dec, respectively. The positional offsets between 
the OGLE and IRSF coordinates were already cor- 
rected. 

IRAF (Image Reduction and Analysis Facility) 5 
software package was used to perform dark and 
flat-field corrections, followed by sky background 
estimation and subtraction. Photometry, includ- 
ing point spread function (PSF) fitting, was car- 
ried out with the DAOPHOT package (Stetson 
1987). We used the DAOFIND task to identify 



4 In paper I, we have observed fields closer to the GC where 
the range of extinction is large (about 3.5 mag in the J 
band), and thus we were able to determine the extinction 
ratio in the J, H, and Kg bands accurately. However, the 
fields observed in this study are located at b ~ —2°, and 
the range of extinction is 0.5 mag in the J band (see the 
bottom panel in Fig. 4), which is less than 0.2 mag in the 
Ks band. 

5 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 



point sources, and the sources were then input 
for PSF-fitting photometry to the ALLSTAR task. 
About 20 sources were used to construct the PSF 
for each image. 

Each image was calibrated with the standard 
star #9172 (Persson et al. 1998), which was ob- 
served every half an hour. We assumed that 
#9172 is J = 12.48 in the IRSF/SIRIUS system. 
The average of the zero point uncertainties and 
the lOcr limiting magnitude in the J band were 
~ 0.02 mag and 16.8 mag, respectively. 

Astrometric calibration was performed, field by 
field, with reference to the positions of stars in 
the 2MASS point source catalog (Skrutskie et al. 
2006). Only the stars with the photometric error 
of less than 0.05 mag in the 2MASS and our cat- 
alog were used for the calibration. The positional 
difference was finally calculated by using the stars 
of < 0.1 mag photometric error, and we obtained 
the standard deviation of the positional difference 
better than O'.'l. 

2.2. Cross-identification of the SIRIUS 
and the OGLE source 

The stars found with IRSF/SIRIUS and OGLE 
have been cross-identified, field by field, using a 
simple positional correlation. The procedure con- 
sists of two steps. First, identification was per- 
formed with large radius (1'.'5) to evaluate astro- 
metric offset between IRSF and OGLE coordi- 
nates. We found that the offsets in R.A. and Dec. 
are in the range from — ft' 7 to +ft'l, and from 
+ft'l to +ft'4, respectively. The offsets seem to 
depend on the position in the OGLE field. Second, 
the OGLE coordinates were corrected for these off- 
sets, and a search radius of ft' 7 was used for the 
identification. 

We show the histograms of positional differ- 
ences of the finally identified stars in Fig. 2. We 
found ~ 25, 000 matches in each field with an rms 
error in the difference between SIRIUS and OGLE 
coordinates 0'.'13 in R.A. and Dec. The values of 
a obtained by fitting with a Gauss function arc 
ft/088 in R.A. and ft/081 in Dec. 

2.3. Data Analysis 

To measure the ratios of total to selective ex- 
tinction Av/Ey-j and Aj/Ey-j, we selected the 
bulge RC stars, which constitute a compact and 
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well-defined clump in a CMD, and are thus good 
tracer of extinction and reddening. The V versus 
V—J CMD constructed with the IRSF and OGLE 
data is shown in Fig. 3. In this analysis we follow 
the procedure described in Udalski (2003), Sumi 
(2004), and paper I. 




V- J 

Fig. 3. — V versus V — J color magnitude di- 
agram of the OGLE bulge fields. Contours are 
from to 1600, linearly spaced by 80. The solid 
and dashed lines show the direction of the inter- 
stellar extinction with Av/Ev-j = 1.251 (this 
work), and 1.393 (Cardelli, Clayton, & Mathis 
1989, Ry = 3.1), respectively. 

First, we divide each field of SIRIUS into nine 
subficlds of <~ 2.'8 x 2f 8 on the sky. Then we con- 
struct V versus V — J and J versus V — J CMDs 
for each subfield. Second, in the CMDs, we extract 
stars in the region dominated by RC stars, and the 
stars are used to make magnitude histograms (lu- 
minosity function, see also Fig. 2 in paper I). The 
peaks of RC stars are fitted with a Gaussian func- 
tion. Third, the stars in the range fitted above 
are employed to see the distribution of RC stars 
in the color V — J, and the color peaks of RC stars 
are also fitted with a Gaussian function. Since the 
mean V magnitudes of RC stars become too faint 
in highly reddened fields, estimates of the peak 
magnitudes and the colors of RC stars can be un- 
reliable in such fields. To avoid this problem, wc 
do not use the subficlds in which we could not 



find a clear peak of RC stars. Note that this ex- 
clusion was made in the second and third steps, so 
a field employed in the J versus V — J CMD can 
be excluded in the V versus V—J CMD if its V 
distribution does not have a clear Gaussian peak, 
and vice versa. As a result, 48 (V versus V — J) 
and 8 (J versus V — J) out of 288 (32 x 9) sub- 
fields were excluded from our analysis. Therefore 
the data points of the two CMD are not identical, 
and the resultant slopes do not necessarily satisfy 
the relation A v /E V -j = Aj/Ey-j + l 6 . 

3. Results 

Combining the IRSF J band and the OGLE 

V band data sets, we have made the V versus 

V — J (top panel in Fig. 4), and J versus V — J 
(bottom panel) CMDs in which the location of 
the RC magnitude and color peaks are shown. 
Error bars include uncertainties in the RC peak 
and photometric calibration. The least-squares 
fits to the data points provide us with the slope 
in the CMDs, Ay/Ey-j = 1.255 ± 0.004 and 
Aj/E V -j = 0.225 ± 0.005. 

The errors in the slopes obtained by using the 
least-squares fit seem to be underestimated be- 
cause the dispersion of the data points is large 
compared with the error bars, which is particu- 
larly noticeable in the J versus V — J CMD. Hence 
we estimated the errors of the slopes by fixing the 
X 2 of the fit to equal the number of degrees of free- 
dom, x 2 /dof = 1, under the assumption that the 
errors are all equal. Application of this procedure 
to the data in both CMDs results in Ay /Ey-j = 
1.251 ± 0.014 and Aj/E v _j = 0.225 ± 0.007. 

These ratios of total to selective extinction pro- 
vide the ratio of total extinction Aj/Ay, with a 
simple algebra of ((A v /E v _j) - 1)/(A V /E v _j) 
and (Aj/E v _j)/((Aj/Ey_j)+l), yielding Aj /Ay = 
0.201 ± 0.011 (from A v /E v _j) and Aj/Ay = 
0.184 ± 0.006 (from Aj/Ey-j). Therefore we 
obtain the weighted mean and error of them, 
Aj/Ay = 0.188 ± 0.005. However, the differ- 
ence of the two values is 0.017, and the combined 
sigma is V0.011 2 + 0.006 2 = 0.013. Hence AjjAy 
could have an error of an order of 0.01. 



'We made a J versus V — J CMD excluding 48 subficlds, 
which are not used in V versus V — J, and confirmed that 
the exclusion of them changes the slope in the CMD by 
only 0.002. 
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Fig. 4. — Location of RC peaks in V versus V — J 
(top) and J versus V — J (bottom) CMDs. The 
solid lines are the least-squares fits to the data 
points. 

4. Discussion 

4.1. Reliability of Our Data Sets 

The faintest data points in the CMDs (Fig. 4) 
are V = 20.2 and J = 14.7. Udalski (2003) set a 
safety margin and regarded his points as complete 
to the limit of V = 20.7, which is still fainter than 
our faintest data point by 0.5 mag. As described 
in section 2, the limiting magnitude in the J band 
is 16.8, also much fainter than J = 14.7. Therefore 
these margins seem to be large enough to have the 
mean magnitude of RC stars reliably measured. 

The peak magnitude of RC stars could be al- 
tered by the dropping completeness. Therefore we 
have checked the completeness and the change of 
the peak magnitudes. First, we have made exper- 



iments in which we added artificial stars of vari- 
ous known magnitudes to our original images, and 
subjected them to the same procedure described 
in §2.1. The detection rates drop as the magni- 
tude becomes fainter; w 95% at J = 14, w 85% at 
J = 15, w 70% at J = 16. Second, we have recon- 
structed luminosity functions with the detection 
rates compensated, and fitted the RC peaks again. 
Comparing the RC peak magnitudes with recon- 
structed ones, we have obtained that the mean 
difference between them is 0.017 mag, and its rms 
is 0.006. The rms is very small compared to the 
error of each data points in the CMDs, typically 
0.02-0.03 mag. We cannot find clear dependence 
of the difference on the peak magnitude, and on 
the number of stars. Hence we can conclude that 
the completeness effect does not change the slopes 
in the CMDs. 

The zero-point uncertainty should be checked 
to derive the slope in CMDs reliably. We ex- 
amined the internal consistency of the duplicate 
sources in overlapping regions of adjacent fields. 
Histograms of mean magnitude difference with the 
sources of photometric error less than 0.05 mag for 
each field set are shown by the left panel in Fig. 
5. The differences of adjacent fields in the direc- 
tion of R.A. and Dec. are shown by white and 
hatched histograms, respectively. We determined 
the rms of the magnitude difference to be less than 
0.02 mag. For another check, we have made com- 
parisons with the 2MASS catalog in the J band. 
The histogram of the mean difference between the 
2MASS and IRSF J magnitudes for each field is 
shown by the right panel in Fig. 5. The mean and 
rms variance of the histogram are 0.009 mag and 
0.015 mag, respectively. The rms variance is sim- 
ilar to the zero-point uncertainty we derived from 
the observation of the standard stars. We there- 
fore conclude that the systematic error in our data 
sets is about 0.02 mag. Note that, in this study, a 
systematic error does not come from absolute zero- 
point uncertainty, but from relative one. Even if 
the absolute magnitudes of RC stars has a sys- 
tematic offset, the resultant plots shown in Fig. 4 
move as a whole with the slope unchanged. Also, 
the right panel in Fig. 5 shows that the absolute 
magnitude offset in the J band is small. 
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4.2. Error Estimates of R\ 
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Fig. 5. — (Top) Histograms of mean magnitude 
differences of adjacent fields with the sources of 
photometric error less than 0.05 mag. The dif- 
ferences in the direction of R.A. and Dec. are 
shown by white and hatched histograms, respec- 
tively. The mean difference and rms of the his- 
tograms are -0.002 and 0.016 (R.A.), and 0.002 
and 0.019 (Dec). (Bottom) Histogram of mean 
magnitude difference between 2MASS and IRSF 
sources in the J band. The mean and rms are 
0.009 and 0.015, respectively. 



The difference of distance to RC stars in the 
observed area could be a factor of the systematic 
error, which become large in the situation that the 
distance is correlated with interstellar extinction. 
However, the distance changes along the Galactic 
longitude due to the presence of the bar structure 
(e.g., Nakada et al. 1991), and the extinction seems 
to change along the Galactic latitude in our ob- 
served area (Sumi 2004), suggesting uncorrelation 
between them. The extent of the observed area in 
the Galactic longitude is only 1?1, which leads to 
an only ~ 0.03 mag difference in the brightness of 
RC stars (Nishiyama et al. 2005), also suggesting 
the small systematic error by the difference of the 
distance. 

Although population effect of the RC stars is 
a source of the systematic error, the dependence 
of the RC brightness on the metallicity is very 
weak, and the metallicity gradient is expected to 
be very small in our small observed area (Udalski 
et al. 2002; Nishiyama et al. 2005). In addition, 
the metallicity should be again correlated with the 
interstellar extinction, and this situation seems to 
be unlikely. 

To examine the systematic errors for the popu- 
lation effect and the different distance of RC stars, 
we made linear fits to the data points in the J ver- 
sus V - J CMD for four OGLE fields separately, 
keeping the same slope (0.225) but changing the 
intersect. The intersects we obtained are 13.37, 
13.37, 13.35, and 13.38 for SC37, 3, SC4, and 
SC39, respectively. Next, to check the dependence 
of the intersect's deviation on the slope value, we 
also made linear fits for the four fields, with chang- 
ing the slope values between 0.211 (0.255— 2a) and 
0.239 (0.255 + 2a). We could not find any system- 
atic trend, and they are within the error of the 
intersect, 13.37 ± 0.03 for all subfields. We thus 
conclude that the difference of the distance and 
population effect do not affect our results. 

The change of the effective wavelength in the 
J band is very small. As shown in the bottom 
panel of Fig. 4, the peak magnitudes of RC stars 
are distributed between mj = 14.1 and 14.7. Us- 
ing the distance modulus of 14.38 toward the GC 
(Nishiyama et al. 2006b) and the absolute magni- 
tude Mj w -0.3 of RC stars (Bonatto et al. 2004), 
Aj can be estimated to be 0.02 < Aj < 0.68. In 
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this range of extinction, the effective wavelength 
in the J band for typical bulge RC stars changes 
by only 0.002^m (see Nishiyama et al. 2006b, §5.1, 
for more details) , and the change of effective wave- 
length is thus negligible. 

4.3. Extinction A\/Ay toward Galactic 
Bulge 

In the wavelength range 0.5 — 0.9/im, extinc- 
tion toward the Galactic bulge has been charac- 
terized by a "steep" curve whose ratio of total 
to selective extinction is small. The steep extinc- 
tion curve was introduced to explain the anoma- 
lous (V — I)o color of RC stars in Baade's window 
(e.g., Popowski 2000; Gould et al. 2001). The 
lower value of the ratio of total to selective ex- 
tinction was also reported from the MACHO V 
and R photometry (Ay / Ey-R ~ 3.5, Popowski et 
al. 2003), and from the OGLE V and / photome- 
try (R VI = A v jE v _i ~ 2.0, Udalski 2003; Sumi 
2004) toward the Galactic Bulge. The empirical 
analytic formula of CCM or Fitzpatrick (1999) 
corresponding to the extinction toward the Galac- 
tic bulge has the single parameter Ry of ~ 2, 
which is much lower than the average value 3.1 for 
diffuse regions in the local interstellar medium. 

The CCM formula uses a simple power law 
^-i.6i m wavelength region > 0.9^m and it 
seems independent of Ry- However, since Ay de- 
pends on Ry, the ratios Ay/Ey-j and Aj/Ey^j 
are dependent on Ry, and these ratios derived in 
our work 7 correspond to Ry ~ 1.8, which is also 
very small. The reddening vectors in the V ver- 
sus V — J CMD for the case of Ry =3.1 and our 
results are plotted in Fig. 3. The figure shows a 
clear difference in the reddening direction. 

Small values of Ry are generally considered to 
indicate the prevalence of small dust grains which 
affect the extinction curve in the ultraviolet - op- 
tical wavelengths. Although a substantial num- 
ber of lines of sight with low Ry values are found 
especially at high Galactic latitude (e.g., Larson 
et al. 1996; Larson & Whittet 2005), only a few 
lines of sight exist with Ry < 2.0 (Szomoru & 
Guhathakurta 1999; Larson & Whittet 2005). 

Next, we try to extend the extinction curve 
to 2/xm. Paper I determined the dependence of 



the interstellar extinction in the J, H, and K$ 
bands toward the GC, and thus we can combine it 
with the result obtained in this paper to determine 
A\j Ay in the H, and K$ bands. Here we should 
take into account the variation of the extinction 
law in different lines of sight, because the area ob- 
served in paper I (| I |< 2?0 and 0?5 <| b |< 1?0) 
does not overlap with that of this study. The vari- 
ation of Ak s I 'Eh-Ks was estimated to be as large 
as ~ 7% in the region of 4° x 2° at the GC (pa- 
per I), and hence we adopt this value as a varia- 
tion of Aj/Ay, re sulting AjJAy = 0.188 ± 0.014 
where 0.014 = ^(0.005)2 + (0.188 x 0.07) 2 . By 
using the ratio Aj : Ah : Ak s = 1 : 0.573 ± 
0.009 : 0.331 ±0.004 (paper I), we obtain Aj/Ay : 
Ah I Ay : A Ks /A v = 0.188±0.014 : 0.108±0.008 : 
0.062 ±0.005. 

Table 1: The wavelength dependence of the inter- 
stellar extinction. 



Ay/Ey-j 
Aj/Ey_j 



Aj/Ay 

Ak s IAv 



IRSF 



/dH^ CCM^ 



1.251 ±0.014 1.325 1.393 
0.225 ±0.007 0.325 0.393 



0.188 ±0.014^ 0.245 0.282 
0.108 ±0.008^ 0.142 0.190 
0.062 ± 0.005^ 0.088 0.118 



7 Rv ~ 1-8 was obtained to reproduce our result Aj/Ay 
0.193 by using the CCM formula. 



References. — (1) van dc Hulst (1946); 

(2) Cardelli, Clayton, & Mathis (1989) for the case of R v = 

(3) The error coming from variation of extinction law in 
different lines of sight is included. 

The resultant A\/Ay in the JHKs wavelength 
range is a steeply declining function. As shown in 
Table 1, the CCM curve (for R v = 3.1), which 
is based on Rieke & Lebofsky (1985), decreases 
much more slowly toward the longer wavelengths. 
In particular, the K$ band extinction is slightly 
greater than one tenth of the visual extinction Ay 
in the CCM curve, which contrasts with the van de 
Hulst (1946) curve, where Ak s is slightly less than 
one tenth of Ay. The derived extinction toward 
the Galactic bulge decreases more steeply as the 
wavelength increases. The steep decrease is rather 
striking, but it was already evident in Messineo et 
al. (2005) and paper I, where a steep extinction 
power law index a ~ 2.0, consistent with the po- 



: 3.1. 
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larization law up to <~ 2.5/im (Nagata et al. 1994), 
was proposed for the GC extinction. To confirm 
this, deep optical and infrared observations in the 
V to Ks bands for the same fields with appropri- 
ate extinction would be important. 
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